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PACS 03.65.Yz - Decoherence; open systems; quantum statistical methods 
PACS 05 . 70 . Ln - Nonequilibrium and irreversible thermodynamics 
PACS 32 . 70 . Cs - Oscillator strengths, lifetimes, transition moments 

Abstract -We study the internal dynamics of an elementary quantum system placed close to a 
body held at a temperature different from that of the surrounding radiation. We derive general 
expressions for lifetime and density matrix valid for bodies of arbitrary geometry and dielectric 
permittivity. Out of equilibrium, the thermalization process and steady states become both quali- 
tatively and quantitatively significantly different from the case of radiation at thermal equilibrium. 
For the case of a three-level atom close to a slab of finite thickness, we predict the occurrence of 
population inversion and an efficient cooling mechanism for the quantum system, whose effective 
internal temperature can be driven to values much lower than both involved temperatures. Our 
results show that non-equilibrium configurations provide new promising ways to control the state 
of an atomic system. 
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Introduction. — Thermalization mechanisms in 
quantum systems driven by changes of external param- 
eters offer a great variety of relaxation phenomena, typ- 
ically studied for many-body systems [ljj5]. What hap- 
pens to the internal dynamics of an elementary one-body 
quantum system in presence of an environment driven out 
of thermal equilibrium? This configuration characterizes 
several systems in biology 6 and physics [7-12 and, de- 
spite its simplicity, it may offer a great richness. Systems 
out of thermal equilibrium have been recently subject of 
intense investigations concerning heat transfer |13]j2l] and 
Casimir-Lifshitz interaction [7| 9 , 20—26 . There, the ther- 
malization dynamics is not considered, being the full sys- 
tem assumed fixed in a given configuration. 

The key quantity for the internal atomic dynamics is the 
lifetime. Approaching an atomic system close to a body 
modifies the local electromagnetic (EM) field surrounding 
it, resulting in lifetimes depending on their relative dis- 
tance, on geometrical and optical properties of the body, 
as well as on the temperature of the system. Typically life- 
time has been studied for configurations at thermal equi- 



librium, when the radiation impinging on the body is at 
thermal equilibrium with it 27-29 . Such studies have 



( a ) E-mail: bruno . bellomo@univ-montp2 . f r 



particular relevance in atom-chip experiments 30,31 and 
in all experiments where the near field of a body is probed 

01211111]. 

Here we study an elementary quantum system offering 
a rich non-equilibrium dynamics: a three-level atomic sys- 
tem placed close to a body of arbitrary geometry and di- 
electric permittivity, held at temperature Tm, and em- 
bedded in an additional radiation coming from the walls 
surrounding the entire system and held at temperature 
Tw (see Fig. [I]). When the walls are irregular and far 
enough from the atom-body system, their contribution can 
be described as a black-body radiation in the region of the 
system IT]. What happens to lifetimes, coherences, and 
steady states? How are these quantities modified with re- 
spect to the thermal equilibrium configuration? What is 
the role of geometric and dielectric properties of the body? 
In this letter, we answer such questions, which have both 
fundamental nature and experimental relevance for sys- 
tems which are naturally out of thermal equilibrium, such 
as recent studies of cold atoms close to superconducting 
surfaces 35,36 , or involving the tip of an AFM close to 
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Fig. 1: The temperature of the body, Tm, and the one of the sur- 
rounding walls located at large distance, T\y, may be different, and 
are kept fixed in time realizing a stationary configuration out of ther- 
mal equilibrium. For the three-level atom, having position R = (r, z) 
and states of increasing energy |1), |2) and |3), we consider a A con- 
figuration with only two allowed dipole transitions between |1) and 
|3) and between |2) and |3), respectively separated by energies ftw^i 
and hui32- 



cold samples 37,38 . Our analysis, concerning both real 
and artificial atomic systems, can be readily extended to 
more complicated level schemes. 

Atomic dynamics. — The atom -field interaction 
is described using the multipolar-coupling Hamiltonian 
Hj = — D-E(R), where D is the atomic electric dipole op- 
erator (supposed to be purely non-diagonal and with non- 
zero matrix elements dis = (1|D|3) and CI23 = (2|D|3)) 
and E(R) is the total electromagnetic field at the atomic 
position. The state of the atom is represented by its den- 
sity matrix whose exact time evolution is governed, 
in the interaction picture, by the trace on the degrees of 
freedom of the field of the von Neumann equation for the 
total density matrix ptot(^) = — j[[Hi(t),pt t(t)]. We treat 
here the atomic dynamics in the limit when Born, Marko- 
vian and rotating- wave approximations hold 1 39 , deriving 
a master equation in terms of transition rates associated 
to the allowed transitions. The validity of these approxi- 
mations is guaranteed by the weak coupling between the 
atom and the field due to the small value of the electric- 
dipole matrix elements. The master equation reads 



d 
dt 



pit) = -i jy^ n |n) (n\ +y^S(-u nrri )\m)(7 

n m,n 

+ ^2S(uj nrn )\n)(n\,p(t)^ 
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where T(-uj nrn ) = 7*j(-o; nm )[d mn ]i[d mn ]J, 

r(cJ nm ) = ^,j7ii(^nm)[d mn ]|[dmn]j, Tl = {1,2,3}, 

(nm) e {(31), (32)} and u 31 ^ cj 32 . S(u nm ) and 
S(—uj nm ) are frequency shifts not playing any role in the 



population dynamics. In the previous equation, Jij(w) is 
defined by 



i 



dse iujs (E l (K,s)E J (Il,0)), (2) 



where we used homogeneity in time for the field correlation 
functions. The steady solution of this master equation 
gives coherences equal to zero and populations 



^>n(oo) N 

P22(00) 



1 

Z 



,(31) 



(1 



,( 31 )^( 32 ) 



(3) 



3n (31 V 32) 

orL ef£ ™eff 



,6 eff 
,(31) 



6 eff 



,(32) 



where, for (nm) e {(32), (31)}, 

(nm) ^(<^nmo 



nm 



with n(uj, T) = (ex-p[hco/k B T] - l) 



(4) 

The coefficients 



n 



(nm) 
eff 



associated to the two allowed transitions satisfy 

(nm) 



^(^nmi -^min ) < n^ ru) < n( 

^nm^max)) where T m 
min{T w ,T M } and T max = max{T w ,T M }. a w (uj nm ) 
and c^M(^nm) are two temperature-independent functions 
strictly connected to the EM field correlators and carrying 
all the dependence on geometrical and material properties 
of the body. The correlation functions of the total field 
appearing in eq. ([2| can be given in terms of the correla- 
tors of the fields emitted by each source by expressing the 
total field in terms of the reflection and transmission op- 
erators 1Z and T associated to the side of body on which 
the atom is located. Each source can be treated inde- 
pendently as if it was at thermal equilibrium at its own 
temperature and the correlators of the source fields can 
be then characterized by using the fluctuation-dissipation 
theorem [21]. The expressions of %(w nm ) and «M(w nm ) 
can be thus obtained, for an arbitrary frequency uo = u nm , 
as a function of 1Z and T as 

3nc^ ^ [d mn ]*[d 
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where i,j G {x,y,z}. A mode of the field is identi- 
fied by the polarization index p =TE,TM, the transverse 
wavevector k = (k x ,k y ), the frequency uj and the direc- 
tion of propagation along the z axis (j) = ±. The three- 
dimensional wavevector is defined as = (k, (j)k z ), its z 

component being k z = \J — k 2 . The polarization vec- 
tors e^ E (k, uj) and e^ M (k, uj) appearing in eq. ([5| are 
defined as e^ E (k,uj) = £(— fc y x + fc^y) and e^ M (k,cj) = 
^(—kz + 0fc^k), where x, y and z are the unit vectors 
along the three axes and k = k/fc. 7^7^ and p^j^ are 
defined as (p, k|^ pw/ew) k') = *£(p, k|n (pw/ew) k'), 
II^ pw ^ and n^ ew ^ being the projectors on the propagative 
(ck < uj, real k z ) and evanescent (ck > uj, purely imagi- 
nary k z ) sectors respectively 21 . 

The steady state of eq. (|3j) gives a thermal state at 
temperature T when, for (nm) G {(32), (31)}, n^ 71 ^ = 
n(unm,T). This condition is trivially verified at thermal 
equilibrium, when Tw = Tm = T. In this case, a peculiar 
cancellation of «wKm) and ^Km) occurs in eq. Q, 
and populations depend only on the ratios fkj nm /&#T, not 
being affected by the body. Remarkably, this is not the 
case if the system is driven out of thermal equilibrium. 
For Tw 7^ Tm the steady state Q does not necessarily 
coincide with a thermal state, and the decay rates can be 
expressed as: 



r(^nm) 
f(-Wnm) 



To(uJ nm ) [«w(w nm ) + a M (Unm)] 
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where r (o; nm ) = ^3™^^ * s tne spontaneous-emission 
rate in vacuum. The quantities T(±uj nrn ) are confined 
between their equilibrium values at T m i n and T max . More- 
over, they equal the ones of thermal equilibrium at an 
effective temperature T^ 771 ^ , associable to each transition, 



defined by n (uj nm , ) 



(nm)x _ (nm) . 



6 eff 



i.e. 



rp(nm) huJ n 

eff ~~k E 



with in general T^ 2 ^ ^ T^\ If the state is not a ther- 
mal one, no particular thermodynamical meaning is asso- 
ciated to the notion of effective temperature. Nonethe- 
less, it is a useful and elegant mathematical tool to de- 
scribe the global dynamics, its meaning being the tem- 
perature for which the decay rates at equilibrium coincide 
with those out of equilibrium. One can then readily inter- 
pret the global dynamics in terms of thermal-equilibrium 
physics, with the strong qualitative and quantitative dif- 
ference that the two transitions feel different temperatures 
whose values depend on the system-body distance, the 
geometry of the body and the interplay of these param- 
eters with the body optical resonances. By varying the 



log 1 



f 1 + (nm) ) 



(7) 



various parameters one can control separately the two ef- 
fective temperatures. One can also obtain T' 



(32) 
eff 



(31) 



eff 



and in this case, even if the full system is out of ther- 
mal equilibrium, the steady atomic state is a thermal one 
at a temperature always between Tw and Tm- This sub- 
tle mechanism allows the emergence of various interesting 
and counterintuitive dynamical features. For example, in- 
version of population ordering of the two lowest-energy 
states |1) and |2) may occur as soon as nf^ < nf^ (i.e. 
as as soon as uj^/T^^ > ^31 /T^ 1 ^). This can happen if 

^(^32,^min) < n(cJ 3 i , T max ) . 

While the quantities r(±uj nrn ), n^ 71 ^ and T e ^ m ^ asso- 
ciated to a given transition (nm) are confined between 
their thermal-equilibrium values at T m i n and T max , this 
is not the case for the steady populations pn(oo) and 
P22(oo). For example, it can be shown that the max- 
imum of pii(oo) (minimum of ^22(00)), obtained when 
n eff } = n (^32,T max ) and n^ 1} = n(uj 31 , T min ), is larger 
(smaller) than its value when Tw = Tm = T min . This 
can be understood since the transition (32) to which the 

(32) 

maximal temperature is associated, T^ ; = T max , is more 
reactive than the transition (31) to which the minimal 



temperature is associated, T e v ff 



(31) 



Equivalent ly, 



the minimum of pn(oo) (maximum of ^22(00)) is obtained 



when T 



(32) 



eff 



T m i n and T eff 



(31) 



Numerical analysis. — In the following we consider 
a specific example providing a numerical investigation of 
the peculiar effects described above for an arbitrary body 
Q As body we consider a slab of finite thickness 5. In this 
case the reflection and transmission operators are diagonal 
in the (k, uj) basis. Their matrix elements are given by the 
Fresnel reflection and transmission coefficients modified by 
the finite thickness 5 of the slab: 



p p (k,uj) = r p (k,uj)- 
r p (k,uj) 



r 2 p (k,uj)e 2ik ^^ 
t p (k,uj)i p (k,uj)e^ kz ^- k ^ 5 



(8) 



r 2 (k,uj)e 



2ik z 



Here r p (k, uj), t p (k,uj) and t p (k,uj) are the ordinary Fres- 
nel coefficients, k zm = \J ^e(uj) — k 2 is the z component 
of the wavevector inside the medium and e(uj) is the di- 
electric permittivity of the body. For this specific case, 
the general formulas for %Km) and au(^nm) in terms 
of reflection and transmission operators reduce to 



i + B(uj nm ) + 2C(uj nm ) ~ 
a-w{uj nm ) = • ct n 

i - B(uj nm ) + 2D(cj nm ) ~ 

OLMK^nm) = ^ ' n 

where d nm = (| [d nm ] x \ 2 , \ [d nm } y \ 2 , | [d nm ] 2 | 2 )/|d, 
i = (1,1,1) and 



(9) 



2 

nm I 1 



1 A11 numerical calculations refer to isotropic dipoles. 
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M+(fc)Im(pp(k,w)), 



(10) 



being M% = (1,1,0) and 



r (^| 2 ,^| 2 ,2fc 2 ). 



In eq. (10), B(cj) and C(cj) are denned in the propaga- 
tive sector while D(cj) in the evanescent one. C(u) and 
D(cj) depend on the atomic position and go to zero for 
large z while, for small z, D(cj) diverges so that «m(^) 
dominates in this region. As a consequence, all the effec- 
tive temperatures tend to Tm in this limit (see eqs. Q 
and ([7])) and thus the atom thermalizes at the body tem- 
perature Tm- B(cj) is independent of z and then, for z 
large enough, it is the dominating contribution in 
and cum(^)- As a result, %(w) > q;m(w) in this region 
and the position of each T e ^ m ^ in the interval [T min , T max ] 
is governed by the value of B(cj). It follows that it always 
exists a distance z for which awKm) = «MKm)- This 
point delimits the two zones of influence where each tem- 
perature dominates for that specific transition. We finally 
observe that at thermal equilibrium B(cj) does not con- 
tribute, since all the quantities are proportional to the sum 
%Km)+«MKm), independent of B(cj) and asymptot- 
ically equal to 1. 

For the numerical simulations, we consider a slab made 
of silicon carbide (SiC). Its dielectric permittivity is de- 
scribed using a Drude-Lorentz model 40 , implying a res- 
onance at uj r = 1.495 x 10 14 rads -1 and a surface phonon- 
polariton resonance at u p = 1.787 x 10 14 rads _1 . A rel- 
evant length scale in this case is c/uo r ~ 2 /am while a 
reference temperature is huj r /kB — 1140 K. 

In Fig. [2| r(-u))/T (uj) is plotted as a function of 
z. The results out of thermal equilibrium are compared 
with the two cases at thermal equilibrium at Tw and 
Tm- We clearly see that at small separation distances 
thermal-equilbrium values at Tm are retrieved, while at 
long separation distances the influence of both Tw and 
Tm is present, as underlined by the arrows indicating the 
asymptotic values. Figure [2] shows that this remains true 
for large values of the slab thickness 5. On the contrary, for 
values of S comparable to the other characteristic lengths 
in the system, the influence of the body is limited with 
respect to z. The asymptotic curves for small z diverge 
as 1/z 3 [28]. The lack of symmetry between the case the 
atom approaches the slab and the case it moves far from 
the body results from the lack of an evanescent term in 
the contribution of the radiation emitted by the walls to 
the local field. The evanescent contribution present in the 
field emitted by the slab is responsible for the fact that for 
z small enough the slab temperature becomes dominant. 



6 = 1 cm 



6 = 3 /im 




0.05 0.1 



Fig. 2: T(—u)/To(lj) as a function of z for uj — 0.5u; r . 
Asymptotic curves for small z (black segments) for Tw = 
T M T max 470K and Tw = Tm = T min 170K. 
Tw = T M = T max (red solid line), (T W ,T M ) (T max ,T min ) 
(purple dot-dashed line), (Tw,Tm) = (T m i n ,T max ) (blue dot- 
ted line), Tw = Tm = T m i n (green dashed line). The symbols 
<(— indicate the asymptotic values (with respect to z) corre- 
sponding to the four couples of temperatures. 



For values of z larger than c/uj ~ 4 /im, the transition rates 
oscillate. This behavior is linked to the integral C(cj) of 
eq. (10), which is purely propagative, and results from the 



phase change of the reflected field 28 . 

As a general remark, we can then state that by increas- 
ing the slab thickness 5 the region of influence of the slab 
temperature becomes larger. The atom-slab distance de- 
termines which temperature will be more relevant in the 
atomic dynamics. This behavior is retrieved in Fig. |3j 
where we plot the effective temperature T^ m ^ as a func- 
tion of z and S for a given frequency co> nm . As previously 



170 200 230 260 290 320 350 380 410 440 470 

10 4 




10 u 

z (|um) 

Fig. 3: Density plot of T^ 171 ^ as a function of z and 5. The 
frequency of the considered atomic transition is uj nm = cj r /2, 
T M = 170 K and T w = 470 K. 
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discussed, for any thickness 5, the atomic temperature for 
small z tends always to the body temperature Tm, while 
for large distances the value of T e ff results from the in- 
terplay between thickness and distance, particularly pro- 
nounced in absence of resonance between atomic transi- 
tion frequency and body optical resonances. Figure [3] also 
shows the existence of regions of oscillatory behavior, orig- 
inating from the propagative part of the spectrum. By 
varying the frequency a different dependence on 6 and z is 
obtained. One can then find couple of frequencies and val- 
ues of S e z such that effective temperatures associated to 
the two transitions are quite different, as explicitly shown 
in Fig. [5] 




Fig. 4: Steady populations as a function of z. For all panels: 0032 = 
ujp and = 2uj r . Left column: 5 = 1 cm. Right column: 5 = 110 
nm. Panels (a)-(b): same color convention of Fig. [2] with T m i n = 
270 K and T max = 540 K. Panels (c)-(d): pn(oo) (red solid line), 
^22(00) (green dashed line) and ^33(00) (dot-dashed purple line) are 
plotted as a function of z in /im for Tw = 540 K and Tm = 270 K. 

In Figs. Qa)-(d) we show that non-equilibrium steady 
populations may exceed their thermal-equilibrium values, 
independent of the atom-slab distance. Moreover, in the 
limit of semi-infinite slab (5 = 1 cm) the populations at 
large distance differ from the ones at thermal equilibrium 
at Tw- In this limit both Tw and Tm contribute to the 
transition rates. On the contrary, at small distances the 
populations always coincide with the ones of a thermal 
state at Tm- The occurrence of ordering inversion of the 
populations pn(oo) and ^22(00) is depicted in panel (d). 
We note that the state |3) is always poorly populated being 
T max smaller than hu p /kB — 1360 K. 

In Fig. [5|a) we highlight the possibility of exploit- 
ing non-equilibrium configurations to provide an efficient 
atomic-cooling mechanism for its internal temperature. 
Indeed, for a large range of atom-slab distances z, the 
steady atomic state is practically a thermal state at a 
temperature lower than the minimal temperature T m i n = 
170 K, as also shown by the inset reporting the distance 
between the atomic state and the closest thermal state, 
Y^Tr(p — a) 2 . In particular, for z = 0.36 /im the atom 
is cooled down to an effective temperature around 48 K 
(temperature of the closest thermal state). It is remark- 
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Fig. 5: For all panels: LJ32 = oo r and UJ31 = uj p . For two values 
of the slab thickness 5, we plot p22(oo) vs pii(oo) by varying z 
(dotted blue line) for T w = 570 K and T M = 170 K. Thermal states 
are represented by solid lines (170 K< T < 570 K: red thick line, 
T < 170 K and T > 570 K: yellow thinner line). In the inset the 
euclidean distance from the closest thermal state is plotted as a 
function of its temperature T c \ by varying z. The distances between 
two thermal states differing of IK, respectively at T = 48 K, 170 K 
and 570 K, are equal to 1.3 x 10 -3 , 1.8 x 10 -3 , 3.4 x 10 -4 . They may 
be useful to estimate how far from thermal states one is by varying 
z. 



able that the above effect can be alternatively produced 
by starting from a thermal-equilibrium configuration with 
Tw = Tm = 170 K, by properly fixing the atomic dis- 
tance and by increasing Tw up to 570 K. As a result, 
the effective atomic internal temperature cools down to- 
wards very low-temperature stationary states. The cool- 
ing mechanism derives from the effective temperatures, 
that at z = 0.36 /im are T^ 2) w 390 K > T e ( ^ 1} w 178 K. 
We remark that the cooling discussed here has nothing to 
do with a reduction of the mean kinetic energy of the atom. 
It refers to the internal atomic temperature associated to 
the Boltzmann factors in the populations. From the in- 
set one sees that at a certain position the steady state 
is exactly a thermal state (the distance from the closest 
thermal state goes to zero). This happens at z ~ 2.25 /im 



where T 



(32) 



eff 



.(31) 
L eff 



520 K. Steady states far from ther- 
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mal ones can also be achieved. This is particularly evident 
in Fig. |5|b) where we show that at z = 0.25 /im the steady 
state is at its largest distance from the curve of thermal 
states with an associated value of p22(oo) of almost 0.9. 
At this position the effective temperatures are inverted, 
with T^ 2) « 227 K and T^ 1] « 476 K. Figure [5] elucidates 
how non-equilibrium configurations provide new tools to 
realize a large variety of different steady states. 

Conclusions. — In conclusion, we have investigated 
the thermalization mechanism of a three-level atomic 
quantum system placed close to an arbitrary body whose 
temperature is different from that of the walls surround- 
ing the atom-body system. We provide closed-form ex- 
pressions for atomic decay rates in terms of the scattering 
matrices of the body valid for arbitrary geometrical and 
material properties in systems both in and out of ther- 
mal equilibrium. We show that configurations out of ther- 
mal equilibrium can be exploited to obtain a large vari- 
ety of steady states, both thermal and non-thermal, with 
populations that can significantly differ from their corre- 
sponding values at thermal equilibrium. Differently from 
the case of thermal equilibrium, they depend on all the 
parameters characterizing the atomic system and its en- 
vironment. In the case the body is a slab, its thickness 
regulates the extension of the zones of influence of the two 
involved temperatures, while the atomic position deter- 
mines which temperature is more relevant in the atomic 
dynamics. Thermalization dynamics can be interpreted 
in terms of effective temperatures associated to each tran- 
sition. We predict peculiar behaviors such as ordering 
inversion of the populations and cooling of the effective 
atomic internal temperature, based on steady configura- 
tions without any additional external laser source. Our 
predictions can be relevant for a wide class of experimen- 
tal configurations involving different physical realizations 
of elementary quantum systems, as real or artificial atoms 
such as quantum dots. A possible experimental realization 
should contain an efficient mechanism to keep the "atom" 
at a given average location. For example, one can imag- 



ine to employ a trapped BEC |9 , or an atomic beam 41 
for ultracold gases, and a mechanical suspension [42] for 
quantum dots. 
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